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ABSTRACT 

The solar cycle and its associated magnetic activity are the main drivers behind changes in the in- 
terplanetary environment and the Earth's upper atmosphere (commonly referred to as space weather). 
These changes have a direct impact on the lifetime of space-based assets and can create hazards to 
astronauts in space. In recent years there has been an effort to develop accurate solar cycle predictions 
(with aims at predicting the long-term evolution of space weather), leading to nearly a hundred widely 
spread predictions for the amplitude of solar cycle 24. A major contributor to the disagreement is the 
lack of direct long-term databases covering different components of the solar magnetic field (toroidal 
vs. poloidal). Here we use sunspot area and polar faculae measurements spanning a full century (as 
our toroidal and poloidal field proxies), to study solar cycle propagation, memory, and prediction. 
Our results substantiate predictions based on the polar magnetic fields, whereas we find sunspot area 
to be uncorrelated to cycle amplitude unless multiplied by area-weighted average tilt. This suggests 
that the joint assimilation of tilt and sunspot area is a better choice (with aims to cycle prediction) 
than sunspot area alone, and adds to the evidence in favor of active region emergence and decay as 
the main mechanism of poloidal field generation (i.e. the Babcock-Leighton mechanism). Finally, by 
looking at the correlation between our poloidal and toroidal proxies across multiple cycles, we find 
solar cycle memory to be limited to only one cycle. 
Subject headings: Sun: dynamo — Sun: activity — Sun: surface magnetism 
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1. INTRODUCTION 

The solar magnetic cycle is without a doubt the main 
driver behind changes in the heliospheric environment 
(Schwenn 2006), violent activity that shapes the Earth's 
magnetosphere (Pulkkincn 2007), and > 80% of the 
Sun's radiative output variability (Domingo et al. 2009) - 
having the potential of disrupting communications, satel- 
lites, and power distribution systems; as well as being 
potentially hazardous to passengers traveling in high- 
altitude polar routes and astronauts in space. Because of 
this, there has been a continuous effort to develop accu- 
rate solar cycle predictions (with aims at predicting the 
long-term evolution of space weather) . 

Cycle predictions are typically classified into three cat- 
egories: extrapolation methods, which use the mathe- 
matical properties of the sunspot data series to predict 
future levels of activity; precursor methods, which use 
different measurable quantities as a proxy to estimate 
the subsequent cycle's amplitude; and model-based pre- 
dictions which use the assimilation of data into models 
of the solar cycle to make predictions (for a review on 
the different types of prediction see for example Petrovay 
2010, and Pesnell 2012). 

As opposed to extrapolation and precursor methods 
(which have been around for nearly 50 years), model- 
based predictions made their first debut at the end of 
solar cycle 23. There were three different model-based 
predictions for cycle 24: two predictions based on mean- 
field kinematic dynamos (Dikpati et al. 2006; Choudhuri 



et al. 2007), and a prediction using a low-mode model 
with variable magnetic helicity which assimilates data us- 
ing an ensemble Kalman filter (Kitiashvili & Kosovichev 
2008). The application of mean-field kinematic dynamos 
to solar cycle predictions has attracted a lot of atten- 
tion because, in spite of using similar models, dynamo- 
based predictions turned out quite different: Dikpati et 
al. (2006) predicted a stronger cycle 24 than cycle 23 
{SSN = 155 - 180), whereas Choudhuri et al. (2007) 
predicted a weaker cycle 24 than cycle 23 {SSN = 80). 
Nevertheless, although similar in nature, the differences 
between these two models are subtle but significant. In 
this work we focus on two of them: 

1. The quantities used to make cycle predictions 
(sunspot area vs. dipolar moment), 

2. The relative importance between diffusive and ad- 
vective transport of the solar magnetic fields (which 
sets the memory of the solar dynamo) . 

One of the main obstacles preventing us from reaching 
consensus regarding these issues is our inability to ob- 
serve the solar magnetic field inside the convection zone - 
leaving surface observations as our only window into the 
inner workings of the dynamo. To compound the prob- 
lem, systematic magnetic field measurements are only 
available since the dawn of the space age (spanning only 
4 cycles) making long-term magnetic proxies necessary. 
The solar cycle can be understood as a process that 
brings the solar magnetic field (back and forth) from a 
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Figure 1. Total daily sunspot area (black dots) is calculated for the northern (top panel) and southern (bottom panel) hemispheres. A 
24-month Gaussian filter is applied to remove the high-frequency component in the data series shown as a solid line for the northern (top 
panel) hemisphere and as a dashed line for the southern (bottom panel) hemisphere. 



configuration that is predominantly poloidal (contained 
inside the r — 9 plane), to one predominantly toroidal 
(wrapped around the axis of rotation; locally perpen- 
dicular to the r — 6 plane). From this follows that at 
least one poloidal and one toroidal field proxies are nec- 
essary. Here we use a recently standardized database of 
polar magnetic flux measurements (based on polar facu- 
lae observations) going back to the beginning of the 20th 
century (Mufioz-Jaramillo et al. 2012), and a long-term 
homogeneous sunspot area database (Balmaceda et al. 
2009), to validate the use of different magnetic proxies 
in the context of cycle prediction (Section S]) , and better 
understand the issue of solar cycle memory (Section [5]) • 
Discussion and conclusions can be found in Section [6l 

2. SUNSPOT AREA: DATA AND SMOOTHING 

In this work we use a homogeneous database of sunspot 
areas (Balmaceda et al. 2009) mainly based on observa- 
tions taken by the Royal Greenwich Observatory, several 
stations belonging to the former USSR (compiled in the 
Solnechniye Danniye bulletin issued by the Pulkovo As- 
tronomical Observatory), and the US Air Force Solar Op- 
tical Observing Network (SOON). We separate the data 
in Northern (top panel in Figure [Ij and Southern (bot- 
tom panel in Figure [T]) hemisphere sets, calculating the 
total hemispheric daily sunspot area. Area belonging to 
groups observed at the equator are not assigned to any 
of the two hemispheres. 

We remove high-frequency components by convolving 
our data series with the modified 24-month Gaussian fil- 
ter: 



Fit,t',a) 



= ^ g-(t-OV(2a^) _ 3~(^*;/)% -2 t'-a<t<t' + a 



t <t' -a 

a<t <t' 
t>t' + a 



(1) 

where t' denotes the position of the center and a = 12 
months the half-width of the Gaussian filter. This type 



of filter has been found to yield more consistent re- 
sults while finding maxima and minima (using different 
databases like the international sunspot number and the 
10.7cm radio flux), than the traditional 13-month run- 
ning mean (Hathaway 2010). 

3. POLAR FACULAE AS A PROXY FOR POLAR 
MAGNETIC FLUX 

Our polar fiux database comes from a recent calibra- 
tion and standardization (Munoz-Jaramillo et al. 2012) 
of four Mount Wilson Observatory data reduction cam- 
paigns (Sheeley 1966, 1976, 1991, 2008). Consecutive 
campaigns were cross-calibrated using five year overlaps 
between different data reduction campaigns (see Fig- 
ure [5]-a) . The cross-calibrated dataset was validated us- 
ing an automatic detection algorithm on intensity data 
from the Michelson Doppler Imager (MDI; Scherrer et al. 
1995) on board the SOlar and Heliospheric Observatory 
(SOHO) spacecraft. The resultant faculae database was 
calibrated in terms of polar magnetic field and fiux using 
magnetic field measurements taken by the Wilcox So- 
lar Observatory and SOHO/MDI (see Figure [2}b). Our 
results demonstrate that there is a strictly proportional 
relationship between facular count, average polar field, 
and total signed polar flux (during all phases of the cy- 
cle and for each hemisphere separately) making it an 
ideal proxy for the evolution of the polar magnetic flelds. 
Once converted into polar fiux values we combine MWO, 
WSO and MDI data into a single consolidated database 
(see Figure [2]-c) and use average polar flux during solar 
minimum as a proxy for the poloidal component of the 
magnetic field (for more details on the validation and 
magnetic calibration of polar faculae data please refer to 
Muiioz-Jaramillo et al. 2012). 

4. RELATIONSHIP BETWEEN THE POLOIDAL AND 
TOROIDAL FIELD PROXIES: IMPLICATIONS FOR 
CURRENT MODEL-BASED PREDICTIONS 

While there is consensus on the process which con- 
verts poloidal to toroidal field (stretching of poloidal 
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Figure 2. Four Mount Wilson Observatory campaigns are standardized using their overlap to obtain a consistent polar faculae database 
(a). Different colors and markers correspond to different data reduction campaigns. This database is calibrated using data from the Wilcox 
Solar Observatory and the MDI, in order into convert it to polar flux estimates (b). The resultant databases are consolidated into a single 
proxy (c). 

flux-transport simulations showing that AR emergence 
and decay leads to the reversal and concentration of po- 
lar flux (Wang, Nash & Sheeley 1989), and which driven 
by semi-synthetic records of sunspot groups obtain a sig- 
nificant correlation between the polar field at minimum 
and the strength of the subsequent sunspot cycle (Jiang 
et al. 2011). Furthermore, a significant correlation has 
been found between AR tilt and the amplitude of the 
next cycle (Dasi-Espuig et al. 2010), suggesting a direct 
connection between AR properties and cycle propaga- 
tion. 

From a practical point of view, the most attractive 
feature of the BL mechanism is the fact that ARs (and 
their role in the evolution of the solar polar magnetic 
field) play a crucial role in the progression of the cycle. 



field by differential rotation; Parker 1955a), there are 
several mechanisms which may be playing a role in the 
conversion from toroidal to poloidal field. Chief among 
them are: the twisting of toroidal field into the poloidal 
plane due to its interaction with helical turbulent convec- 
tion (Parker 1955a), tilted active region (AR) emergence 
and decay (also known as the Babcock-Leighton mech- 
anism, Babcock 1961; Leighton 1969), flux-tube insta- 
bilities (Schmitt 1987), and hydrodynamical shear insta- 
bilities (Dikpati & Oilman 2001) - for a comprehensive 
review please refer to Charbonneau (2010) and references 
therein. 

Currently, the main contending theory for the genera- 
tion of poloidal field is the Babcock-Leighton (BL) mech- 
anism. Supporting evidence comes in the form of surface 
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Figure 3. Correlation between polar flux at solar minimum and 
the amplitude of the next cycle. Square (circular) markers de- 
note data for the northern (souther) hemispheres. Markers are 
numbered using cycle amplitude as reference. The dashed line cor- 
responds to a linear fit using the least absolute residuals method. 
The text inside the figure panel indicates the Pearson's correlation 
coefliicient and its statistical significance. 

allowing modelers to use surface magnetic field observa- 
tions to drive dynamo-based predictions. In particular, 
Dikpati et al. (2006) used sunspot area to drive the gen- 
eration of poloidal field in their model-based prediction, 
whereas Choudhuri et al. (2007) bypassed the generation 
of poloidal field by directly using the axial dipole moment 
component of the solar magnetic field calculated from the 
polar fields by Svalgaard et al. (2005). The axial dipole 
moment at minimum has been found to be a good pre- 
dictor for cycle amplitude; both when calculated using 
direct measurements of the polar fields, for cycles 22-24 
(Schatten 2005; Svalgaard et al. 2005), and when esti- 
mated using solar open magnetic flux derived from the 
historical aa index (for cycles 12-23), after removing the 
contribution of the solar wind speed (Wang & Sheeley 
2009). Taking advantage of our century of toroidal and 
poloidal proxies, our first task is to look at their relation- 
ship from the point of view of cycle prediction. 

We find a good correlation between polar flux at min- 
imum and the amplitude of the next cycle (Figure [3]), 
with a Pearson's correlation coefficient oi p — 0.60 and 
P = 99% confidence level. An interesting feature of the 
correlation between polar fiux at minimum and the next 
cycle's amplitude is the apparent existence of two sepa- 
rate branches in their relationship. This feature becomes 
more evident after performing a linear fit using the least 
absolute residuals method (shown as a dashed line in 
Figure |3]) which naturally gives less weight to possible 
outliers in the dataset. The separation of data into two 
branches results in a large improvement in the Pearson's 
correlation coefficient: p = 0.96 (p = 0.95) and P — 99% 
[P = 99%) confidence level for the main (secondary) 
branches. This suggests that finding a way to evalu- 
ate on which branch will a cycle fall would result in a 
highly effective method of prediction. An in-depth study 
of this separation (which seems to be related to the rel- 
ative strength of the dipolar and quadrupolar moments 
during minimum) and its application to solar cycle pre- 
diction can be found in Muiioz-Jaramillo et al. (2013). 

Looking at the relationship between our toroidal and 
poloidal field proxies we find no correlation between max- 



imum cycle sunspot area (i.e. cycle amplitude) and po- 
lar flux (p = 0.16 & P = 50%; Fig. |la) nor between 
total cycle sunspot area and polar flux (p = 0.19 & 
P = 59%;Fig. ID-b). Although the apparent disconnec- 
tion between sunspot area and polar flux at minimum 
could be interpreted as evidence against the BL mech- 
anism, it is important to highlight that the systematic 
tilt presented by ARs (Hale et al. 1919) is a crucial com- 
ponent of this mechanism of poloidal field generation. 
Indeed, as shown in Figures |4]-c & d, the multiplication 
of maximum cycle sunspot area, and total cycle sunspot 
area, by the area-weighted average tilt (normalized using 
latitude of emergence and calculated from Mount Wil- 
son Observatory data; see Dasi-Espuig et al. 2010) turns 
them into quantities which are correlated to polar flux at 
solar minimum (with p = 0.74 & P = 99% and p = 0.67 
Sz P = 99% respectively). This agrees with the results 
of Kitchatinov & Olemskoy (2011) who find, for cycles 
19-21, a good correlation between the total contribution 
of all sunspot groups to the large scale dipolar field (cal- 
culated using the area of the largest spot, the angular 
extent, and the tilt of each group) and the amplitude of 
the next cycle. 

These results suggests that sunspot area alone may not 
be an appropriate quantity for use in model-based pre- 
dictions and explains partly why the models of Dikpati 
et al. (2006) and Choudhuri et al. (2007) yielded such 
different predictions. However, the improvement intro- 
duced by including tilt suggests that future model-based 
predictions which rely on AR data will likely obtain bet- 
ter performance if they assimilate tilt, as well as sunspot 
area data. It is important to note that these results are 
obtained using a simplified way of including information 
on the spatio-temporal distribution, time and latitude of 
emergence, of ARs (as performed by Dasi-Espuig et al. 
2010). A detailed assimilation of this information (only 
possible using a more sophisticated model) leads to a 
significant improvement in the estimation of solar mini- 
mum conditions (see Cameron et al. 2010 and Cameron 
& Schiissler 2012), and will likely be an integral compo- 
nent of future model-based predictions. 

5. OBSERVATIONAL STUDY OF SOLAR CYCLE MEMORY 

The difference between the model-based predictions of 
Dikpati et al. (2006) and Choudhuri et al. (2007) that 
has received the greatest amount of attention is the rela- 
tive importance of advective and diffusive transport. On 
the one hand, Dikpati et al. (2006) use a low-diffusivity 
model in which the meridional flow - a poleward plasma 
flow observed in near surface layers (Komm et al. 1993; 
) which is believed to turn around into an equatorward 
flow near the bottom of the convection zone, driving the 
equatorial migration of active latitudes (Choudhuri et 
al. 1995) - is the most important mechanism of mag- 
netic flux-transport. On the other hand, Choudhuri et 
al. (2007) use a high-diffusivity model in which diffu- 
sion is the most important mechanism of magnetic flux- 
transport. 

In an influential theoretical study of advection- 
dominated (AD) vs. diffusion-dominated (DD) model- 
based predictions, Yeates et al. (2008) found that the 
main difference between AD and DD predictions is the 
memory-span of the solar dynamo. To reach this conclu- 
sion they performed simulations in which the source of 
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Figure 4. We finti no correlation between the maximum (a), or total (b), sunspot area of a cycle and polar flux during the subsequent 
minimum. However, both maximum (c) and total (d) sunspot area become correlated with polar flux if multiplied by the area-weighted 
average tilt normalized by latitude of emergence (AWNL). Square (circular) markers denote data for the northern (souther) hemispheres. 
Markers are numbered using cycle amplitude as reference. The dashed line corresponds to a linear fit using the least absolute residuals 
method. The text inside figure panels indicate the Pearson's correlation coefficients and their statistical significance. 



poloidal field varied stochastically with time, and looked 
at the correlation between polar flux at the minimum 
of cycle n and the amplitude of cycle n, n + 1, n + 2, 
and n + 3. They found that in the DD regime polar 
flux at minimum is only correlated with the amplitude 
of the following cycle (n + 1), whereas in the AD regime 
polar flux at the minimum of cycle n is correlated with 
the amplitudes of cycle n, n + 1, and n + 2. In a re- 
cent revision of this work including the effect of turbulent 
pumping - magnetic transport associated with the mor- 
phological asymmetry between convective upflows and 
downflows (Tobias et al. 2001) - Karak & Nandy (2012) 
found that the addition of turbulent pumping removes 
any long-term solar cycle memory, turning AD and DD 
predictions undistinguishable. Taking advantage of our 
long-term poloidal and toroidal proxies, we study cycle 
memory from an observational point of view. 

Figure [5] displays the correlation of polar flux at the 
minimum of cycle n and the amplitude of cycles n (a), 
n + 1 (b), n-|-2 (c), and n + 3 (d), showing only significant 
correlation between polar flux at the minimum of cycle n 
and the amplitude of the next cycle (n -I- 1, Fig. [5]-b). In 
the light of the theoretical studies of Yeates et al. (2008) 
and Karak & Nandy (2012), these results suggest that 
purely advection-dominated models of the solar cycle, in 
which meridional flow is more important than turbulent 
mechanisms of flux-transport (i.e. turbulent diffusivity 
and pumping), are inconsistent with observations. How- 
ever, it is possible for random variations in meridional 



flow amplitude to also reduce cycle memory (removing 
this inconsistency). 

6. DISCUSSION AND CONCLUDING REMARKS 

In a recent theoretical study of the dynamo basis 
of precursor predictions, Charbonneau & Barlet (2011) 
found that the polar field has precursor value (i.e. it 
is well correlated with the next cycle's amplitude) only 
when there is a connection between the surface and in- 
terior layers. The correlation we find between polar flux 
at minimum and the amplitude of the next cycle demon- 
strates the need for such a connection between the sur- 
face magnetic field and the bottom of the convection 
zone. Although high-resolution observations of the polar 
field show concentrated patches of magnetic field strong 
enough to be buoyant (Shiota et al. 2012), instead of a 
diffuse large-scale unipolar field, our results suggest that 
the total polar magnetic fiux can still be seen as the sur- 
face manifestation of low-order moments in the multi- 
pole expansion of the solar magnetic field. Perhaps, 
rather than direct subduction of polar magnetic fields, 
fiux-transport mechanisms act on the roots of polar fiux- 
tubes, pinning and stretching them across the bottom of 
the convection zone to form a large-scale poloidal field 
from which toroidal field is inducted. 

The correlations we find between the toroidal and 
poloidal proxies (Figures |3] & Hl-c & d) represent strong 
observational evidence in favor of the BL mechanism as 
the main source of poloidal field in the Sun - by link- 
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Figure 5. Correlation between polar flux at the minimum of cycle n and tfie amplitude of cycle n (a), n + 1 (b), n + 2 (c), and n + 3 
(d). Polar flux at minimum is only correlated to the amplitude of the next cycle (n + 1; b). This set of correlations is only consistent 
with a solar dynamo in which turbulent effects (diffusion and/or pumping) are as important as flux-transport by meridional flow. Square 
(circular) markers denote data for the northern (souther) hemispheres. Markers are numbered using cycle amplitude as reference. The text 
inside figure panels indicate the Pearson's correlation coefficients and their statistical significance. 



ing together the determination of polar fields by active 
region emergence and decay, and its subsequent shear- 
ing by differential rotation to set the amplitude of the 
next cycle. Together they support the current solar cy- 
cle logic where the two components of the solar magnetic 
field (toroidal and poloidal) generate each other sequen- 
tially (nicely illustrated in Charbonneau 20f0). This has 
important consequences from the point of view of cycle 
prediction since it links the entire cycle causality to ob- 
servable quantities in the photosphere. In particular, our 
results substantiate proposed precursor methods based 
on the polar magnetic field at solar minimum (by extend- 
ing their applicability to the last century) and suggests 
that they will be the most successful predictions once we 
pass solar maximum. 

Taking advantage of a century of poloidal and toroidal 
field proxies we analyzed two of the main differences be- 
tween the dynamo-based predictions of Dikpati et al. 
(2006) and Choudhuri et al. (2007) from an observa- 
tional point of view. In the case of data assimilation 
we find that sunspot area is not well correlated with the 
amplitude of the next cycle, unless it is multiplied by 
area- weighted average tilt (see Fig. |4]). This suggests 
that tilt, which is a crucial component of the Babcock- 
Leighton mechanism, needs to be used (together with 
sunspot area or active region flux) in order to obtain 
more accurate model-based predictions. 

In the case of solar cycle memory we find observations 
to be consistent with a short term memory (also found 



by Solanki et al. 2002). This result is only naturally con- 
sistent with diffusion-dominated dynamos, or dynamos 
with significant transport by turbulent pumping (as pro- 
posed by Guerrero & de Gouveia Dal Pino, 2008). How- 
ever, a purely advection-dominated dynamo driven by a 
randomly fluctuating meridional flow may be able to re- 
duce cycle memory in accordance with observations. A 
detailed study of cycle memory in advection-dominated 
dynamos is necessary in order to address this discrep- 
ancy. 
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